What Determines the Yield Stress in Amorphous Solids? 
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A crucially important material parameter for all amorphous solids is the yield stress, which is the 
value of the stress for which the material yields to plastic flow when it is strained quasi-statically 
at zero temperature. It is difficult in laboratory experiments to determine what parameters of the 
inter-particle potential effect the value of the yield stress. Here we use the versatility of numerical 
simulations to study the dependence of the yield stress on the parameters of the inter-particle 
potential. We find a very simple dependence on the fundamental scales which characterize the 
repulsive and attractive parts of the potential respectively, and offer a scaling theory that collapses 
the data for widely different potentials and in different space dimensions. 



Introduction At small external strain every solid re- 
acts elastically. We are interested here in the response 
of amorphous solids to high external strains, and in par- 
ticular in the yield stress cry which cannot be exceeded 
without effecting a plastic response which typically leads 
to mechanical failure via plastic flow, shear banding or 
fracture [IH3]- I n this Letter we focus on the funda- 
mental microscopic features which determine the yield 
stress, using to great advantage the versatility of nu- 
merical simulations in which the inter-particles potential 
can be varied at will. We work at zero temperature and 
quasi-static external straining conditions, (the so-called 
athermal quasi-static or AQS limit) where very precise 
simulation results can be obtained. To introduce the is- 
sue examine Fig. [T] which exhibits a typical stress vs. 
strain curve, here for a system of N = 500 particles in 
AQS conditions. The yield stress is the average stress in 
the elasto-plastic steady state that is obtained, say, after 
100% deformation. 

Clearly, the connection of the yield stress to the fun- 
damental microscopic characteristics of the material is 
of great interest, but in fact far from being determined. 
To a large extent the determination of the yield stress of 
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FIG. 1: Color online: a typical stress vs strain curve. At 
low strains the stress increases linearly with the slope being 
the shear modulus. At a strain value of the order of 10% the 
system undergoes a yielding transition [4(. ay is the mean 
stress computed at the elasto-plastic steady state obtained 
after, say, 100% strain (see red curve). 



various amorphous solids is still a question of trial and 
error without much theoretical input. From the point of 
view of fundamental theory we have results for the depen- 
dence of the yield stress on the density p only in the case 
of very simple inter-particle potentials which are purely 
repulsive, without an attractive tail. Simple results are 



available for potentials of the form 
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ing then the re-scaled (dimensionless) density as p 
(where A is a typical length scale, cf. Eq. J2])), then the 
yield stress for solids at different dimensionless densities 
p was shown to vary like 
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It had been shown [5[ that this scaling law is not valid 
in general when the inter-particle potential contains an 
attractive part, except at very high densities when the 
particles are squeezed against the repulsive part of the 
potential. The aim of this Letter is to provide a scaling 
theory that remains valid for generic potentials with an 
attractive and a repulsive part. 

System and potentials: We investigate a binary 
system where the amount of bi-dispersity of 'small' and 
'large' particles was chosen to guarantee that the models 
produce good glass formers both in 2 and 3 dimensions. 
We constructed a template for inter-particle potentials 
where the microscopic lengths can easily be tuned: 
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where r m { n /X i j is the length where the potential attain 
it's minimum, and r co /Xij is the cut-off length for which 
the potential vanishes. The coefficients a, b and c-it are 
chosen such that the repulsive and attractive parts of 
the potential are continuous with two derivatives at the 
potential minimum and the potential goes to zero con- 
tinuously at r co / Xij with two continuous derivatives as 
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FIG. 2: Color online: a selection of different potentials used to 
determine the yield stress. Panel A depicts potentials where 
rco/Aij is the changing length. Panel B shows potentials for 
which r co /Xij remains unchanged but have a different attrac- 
tive behavior. 
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FIG. 3: Color online: stress vs. strain curves in 2-dimensional 
systems with the potential shown in Fig. [2l panel A. Note the 
convergence with increasing r co . 



well. To satisfy the latter constraints it suffices for n to 
be equal to 3. In other cases where r co is fixed and we 
want to control the shape of the attractive part of the po- 
tential, n = 4. The interaction length-scale Ay between 
any two particles i and j is Ay = 1.0A, Ay = 1.18A and 
Ay = 1.4A for two 'small' particles, one 'large' and one 
'small' particle and two 'large' particle respectively. The 
unit of length A is set to be the interaction length scale 
of two small particles, e is the unit of energy and ks = 1. 

Methods: The work presented here investigates sys- 
tems under simple shear in the athermal quasi-static 
(AQS) limit, T — > and 7 — s> 0, where 7 is the strain 
rate. AQS methods have been extensively used recently 
(fj|-[ll| as a tool for investigating plasticity in amorphous 
systems. The order in which the limits T — ¥ 0, 7 — > 
are taken is important, since one expects that at any fi- 
nite temperature the stress in the system can thermally 
relax given long enough time [l2j (or small enough strain 
rates), hence the limit T — s> should be taken prior to 
the 7 — > limit. According to AQS methods, starting 
from a completely quenched configuration of the system, 
we apply an affine simple shear transformation to each 
particle i in our shear cell, according to 
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in addition to imposing Lees-Edwards boundary condi- 
tions [I3j . The strain increment £7 plays a role analogous 
to the integration step in standard MD simulations. We 
choose for the discussed systems £7 = 5 x 1CP 5 , which 
is sufficiently small for the analysis of the steady state 
mean values. The affine transformation Eq. ([3]) is then 



followed by the minimization |14| of the potential energy 
under the constraints imposed by the strain increment 
and the periodic boundary conditions. 

The yield stress in 2-dimensions: We pre- 
pared systems with density p ranging in the interval 
[0.775, 0.935] with increments of 0.02. The low bound- 
ary of this interval is determined by retaining positive 
pressure in our simulation box for any of the considered 
microscopic lengths. For lower densities the system be- 
comes porous with patches of vacuum, and we do not 
investigate such states. The range of the microscopic 
lengths considered to generate the various potentials is 



1.2Ay x 1.06" 



where k G {0, 12}. (4) 



The microscopic lengths used to generate our data-set 
of cry values include all combinations of r co and p with 
n = 3. A few examples of the different potentials can be 
seen in Fig. [2] panel A. Each such system with N = 500 
particles was strained in AQS conditions and a typical 
result for stress-strain curves with given p is shown in 
Fig. [3] The yield stress in simulations are read from the 
steady-state stress-strain curve as depicted in Fig. [TJ To 
understand these ay values, we note that for a finite cut- 
off length the density determines how many neighboring 
particles are within the interaction length, and ay is ex- 
pected to be a strong function of the density. In the limit 
r co — > 00 all the particles are within the interaction range 
and the dependence on the density disappears in favor of 
a limiting value of ay — > ay . 

To collapse the data on a single graph we note that 
there are three independent lengths in this problem, 
namely A, r co and p~ x l d . These three lengths render two 
ratios, i.e. s\ = X/p^ 1 ^ = p 1 ^ and S2 = r co /p In 
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FIG. 4: Color online: the scaling function /(s) (cf. Eq. [6|, 
obtained by collapsing the computed ay values in 2-D. Hot 
colors are for shorter values of r co , small symbols represent 
lower densities. Potentials used here are of the type found in 
Fig. OH panel A. 

general one expects a scaling law with these two variables 

s 

°y = -raP v 9(si,s 2 ) ■ (5) 

We expect that the limit s\ — ¥ oo should reproduce the 
high density limit Eq. ([T}. In contrast, we expect that 
S2 — > oo should lead to the converged large interaction 
limit ay. Thus these two dimensionless numbers work 
in opposite direction, and we consider the dimensionless 
product s = S] x S2 in which the effects of both ratios 
can balance. These consideration lead to a scaling ansatz 
in the form 

Here the scaling function /(s) is expected to behave 
according to 

p* ( const for s — * oo 
/(S) = ( a-" fors^O. ' (7) 

Indeed, plotting ay / 1 p v as a function of s leads to a su- 
perb data collapse, as one can see in Fig. [4] 

It should be noted that the choice of the scaling vari- 
able s = s\ x S2 can be directly validated by checking 
that the value of ay / p v is indeed invariant for a given s 
for any si and S2 = s/si, see Fig. [SJ In the figure we plot 
the scaling function f(s — s± x s 2 ) vs. f(s = s± x S2) for 
a series of different values of s such that s\ x s 2 ~ si X S2- 

The yield stress in 3-dimensions: A sensitive test 
of the above mentioned scaling theory is provided by 
validating the explicit dimensionality dependence. To 
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FIG. 5: Color online: plotted are the values of f(s) vs. f(s) 
from the dataset presented in Fig. 21 where - w 1. 




FIG. 6: Color online: the scaling function f(s) (cf. Eq. [6|, ob- 
tained by collapsing the computed ay values in 3-dimensions 
using potentials as in Fig. [2]panel A and N=1000. Hot colors 
are for shorter values of r co , small symbols represent lower 
densities. 

this end similar data was collected for samples in 3- 
dimensions using the potentials with n = 3 and vary- 
ing r co values, see Eq. combined with densities in 
the interval [0.76, 1] with increments of 0.03 (emanating 
from the above-mentioned considerations) . The beautiful 
scaling collapse in Fig. [BJ simply changing the numerical 
value of v as dictated by Eq. ([1]) confirms the choice of 
scaling function, Eq. ([6]). 

The effect of potential shape on the yield stress 
in 2-dimensions: We have so far been successful in 
identifying the crucial microscopic lengths involved in 
determining the yield stress. The role of a change of 
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shape in the attractive part of the potential now remains 
to be elucidated. (However, regardless of the shape of 
the potential we expect the limiting values to remain un- 
changed. In other words, we still expect that the limit 
si — > oo should reproduce the high density limit and that 
S2 — > oo should lead to the converged large interaction 
limit 5y.) To this end we used two variations of poten- 
tials with n = 4 and r co as in Eq. (j4]). For examples of 
such potentials see Fig. [2j panel B. 

The change in physics as a result of changing the shape 
of the attractive part of the potential can be easily un- 
derstood. As the potential well widens, the slope of the 
potential on the way to r co increases. As long as the sec- 
ond shell of particles is outside the interaction range, the 
system softens as a result of the widening well. As soon 
as r co reaches a value that includes the second shell, the 
strong attractive forces stiffen the system, thus creating 
a dip in the scaling function. 

To corroborate this picture we repeated the AQS 
straining experiments for the altered potentials (see lower 
panel of Fig. [2] at similar densities as before in 2- 
dimensions. In Fig. [7]wc sec the unchanging asymptotes 
of the scaling function and the expected dip in ay j p v . 
Note the larger dip in the right panel of Fig. [7] Note 
that this effect of softening and then hardening is always 
present even if it is barely noticeable 
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FIG. 7: Color online: the scaling function f(s) (cf. Eq. [6|, 
obtained by collapsing the computed ay values in 2-D. Hot 
colors are for shorter values of r co , small symbols represent 
lower densities. Left panel is the collapse of ay values com- 
puted using the 'dash-dot' fashion potentials appearing in Fig. 
[2]panel B and the right panel is the collapse of ay values com- 
puted using the 'dashed' fashion potentials in the same panel. 

Concluding remarks: We have presented a scaling 
theory of the yield stress in AQS conditions, stressing 
the data collapse, and therefore of the high predictabil- 



ity of the scaling theory, when the potentials are changing 
widely both in their shape and in the cutoff scale. The 
data collapse means that measuring the data for one po- 
tential we can predict how the yield stress will change for 
other potentials. Together with the recent scaling the- 
ory for the flow stress as a function of temperature and 
strain rate (albeit at this point only for purely repulsive 
potentials [5J]), we begin to see a powerful theory based 
on scaling concepts that is emerging for the discussion of 
plasticity in amorphous solids. 
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Note that slight increase in ay for the largest value of 
r co in Fig. This is due to the influence of second shell 
neighbors entering the interaction range. 



